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© Different kinds of video noise filters are known. Normally such filters need a certain amount of memory, e.g. 
field memories, or are not suited for every noise level or are not adaptive to the local structure of the picture. 
The inventive solution is globally and locally adaptive. A noisy input signal x is filtered with a restoration filter of 
median type to generate a filtered input signal y. The sum of the absolute differences between filtered and 
unfiltered signal is calculated for each position of a sliding window within the input signal, representing a local 
estimate of the noise, and is combined with a global measure of the input signal noise to compute two 
coefficients a and b which are respectively applied to the unfiltered and filtered signal to generate the output 
signal z = a"x + b"y. 

Advantageously different kind of filters operate in parallel, whereby the kind of filter elected is locally adapted to 
the picture activity. 
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The present invention relates to a method and to an apparatus for noise reduction. 
Background 

s Different kinds of video noise filters are known. M. Unser, "Improved restoration of noisy images by 

adaptive least-squares post-filtering", Signal Processing 20 (1990), pp. 3-14, gives an example. Normally 
such filters need a certain amount of memory, e.g. field memories, or are not suited for every noise level or 
are not adaptive to the local structure of the picture. 

10 Invention 

It is one object of the invention to disclose a method for noise reduction, which does need only simple 
hardware and takes into account different noise levels and picture structures. This object is reached by the 
method disclosed in claim 1 . 
is It is a further object of the invention to disclose an apparatus which utilises the inventive method. This 

object is reached by the apparatus disclosed in claim 9. 

For the inventive noise reduction the main requirements listed below are considered: 
" field processing i.e. without field memory; 
" few memories and simple hardware realisation; 
20 " processing adaptive to any level of noise; 

" processing adaptive to the local structure of the picture, 
i.e. detecting plain areas, edges and fine details so as to preserve the relevant information of the picture. 
Different filters could be used: 
1 . linear low pass filter with different coefficients; 
25 2. median filter; 

3. linear filter in the direction of the edges (directional filter); 

4. directional median filter; 

5. weighted averages of filters 1, 2, 3 and 4; 

6. filtering with one of the filters 1*. 2, 3 and 4 with fixed threshold and local measure of picture activity by 
30 the Laplacian. 

These filter types have the following properties; 

* median filters are more suitable to remove noise in plain areas than linear filters especially in the case 
of high frequency noise of high amplitude; 

* median filters preserve edges and linear filters do not; 

35 * directional filters are very efficient for preserving edges and picture accuracy while filtering slightly but 
efficiently enough; 

median filters completely lose the fine textures (of a zone plate test picture, for example) and the only 
way to preserve them is to use a directional filter; 

weighted averages of filters give good results but the weighting coefficients are not easy to adapt 
40 locally without correct noise estimation and can lead to hard switching and to the introduction of 

artefacts; 

thresholds and local measures of activity give good results when adapted to the picture but they have 
to be combined with noise estimation on the whole picture. 
The inventive solution is globally and locally adaptive. The implemented filter is derived from that 
45 proposed by Unser (cited above) and works as follows: 

The noisy input signal x is filtered with a restoration filter of median type to generate a filtered input signal 
y. The sum of the absolute differences between filtered and unfiltered signal is calculated for each position 
of a sliding window within the input signal, representing a local estimate of the noise, and is combined with 
a global measure of the input signal noise to compute by least squares regression two coefficients a and b 
so which are respectively applied to the unfiltered and filtered signal to generate the output signal z = a*x + 

b-y. 

Advantageously the only required prior information is the noise variance a 2 , assuming that the original 
signal u is degraded by additive stationary noise n. The size of the local noise estimation window can be 
different from the size of the window used by the restoration filter, wherefore they can be adjusted 
55 independently for optimal performance. 

The weighting coefficients a and b are made optimal in the sense that they are minimising the quadratic 
error t 2 between the output signal z and the noise free original signal u. It is true that the original signal is 
not known at a receiver, but global noise statistics can be calculated at the receiver side (e.g. described in 
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EP-A-92400785) and used for that optimisation. Preferably a + b = 1, so that no bias is introduced by the 
filtering. This solution usually performs slightly worse, but is rather simple to implement, no computational 
complexity need be added. 

If the filtering tends to degrade the signal, a predominant weight will be given to the unfiltered signal x. 
5 Conversely, when the sum of differences is small and close to a reference value corresponding to the 
residue of the noise alone, the weight will be shifted to the filtered signal. 

The choice of the kind of restoration filter does heavily influence the merit of the inventive noise 
reduction. Therefore, advantageously different kind of filters can operate in parallel, whereby the kind of 
filter elected has to be locally adapted to the picture activity. That branch which results in a minimum error 
w between the filtered and the unfiltered signal will be selected. 

In principle the inventive method is suited for noise reduction, in which for estimating a local noise value 
within a sliding window of a picture the quadratic error between an unfiltered noisy input signal and the 
filtered input signal is computed and used for calculating a weighted average, of the filtered and unfiltered 
input signal as output signal, whereby for calculating said weighted average also a global noise value is 
75 used and for computing said local noise value pixel difference signals of said window are calculated. 

Advantageous additional embodiments of the inventive method are resulting from the respective 
dependent claims. 

In principle the inventive apparatus includes: 

- a first filter, a second filter and a third filter having a common input signal and a different kind of 
20 characteristic; 

- a first control circuit controlling the weighting values of a first multiplier operating on said input signal 
and a second multiplier operating on the output signal of said first filter, whereby the outputs of these 
multipliers are combined in first adding means; 

- a second control circuit controlling the weighting values of a third multiplier operating on the output 
25 signal of said first filter and a fourth multiplier operating on the output signal of said second filter. 

whereby the outputs of these multipliers are combined in second adding means: 

- a third control circuit controlling the weighting values of a fifth multiplier operating on the output signal 
of said second filter and a sixth multiplier operating on the output signal of said third filter, whereby 
the outputs of these multipliers are combined in third adding means; 

30 - first subtraction means which feed the difference signal between said input signal and the output 
signal of said first filter to said first control circuit; 

- second subtraction means which feed the difference signal between said input signal and the output 
signal of said second filter to said second control circuit; 

- third subtraction means which feed the difference signal between the output signal of said first filter 
35 and the output signal of said second filter to said second control circuit; 

- fourth subtraction means which feed the difference signal between said input signal and the output 
signal of said third filter to said third control circuit; 

- fifth subtraction means which feed the difference signal between the output signal of said second filter 
and the output signal of said third filter to said third control circuit; 

40 - further subtraction means which feed the respective difference signals between said input signal and 
the output signals of said adding means to a control and mixing circuit which selects an output signal 
of said adding means according to the minimum error derived from these difference signals and 
according to a global noise variance value, 
whereby for each of the branches said weighting values are derived from error signals calculated in the 
45 respective control circuit from said difference signals and from said global noise variance value. 

Advantageous additional embodiments of the inventive apparatus are resulting from the respective 
dependent claim. 



so 



Drawings 



Preferred embodiments of the invention are described with reference to the accompanying drawings, 
which show in: 

Fig. 1 first embodiment of an inventive noise reduction circuit; 
Fig. 2 filter-related directions in the picture; 
55 Fig. 3 window for noise filtering; 

Fig. 4 second embodiment of an inventive noise reduction circuit; 
Fig. 5 correction terms for the comparison of the filters errors. 
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Preferred embodiments 

In Fig. 1 noise n is added in a channel, e.g. a TV channel, to a noise free original signal u., e.g. a TV 
signal. The resulting noisy input signal x of the inventive apparatus, e.g. in a TV receiver or a VCR. is 

s filtered with a restoration filter MF of median type to generate a filtered input signal y. The difference A 
between signal x and signal y is calculated in a subtractor SUB and fed to a control circuit CON which 
forms from A the sum of the absolute differences between filtered and unfiltered input signal for each 
position of a sliding window within the input signal x, representing a local estimate of the noise. 
The global noise statistics can be calculated in a noise measurement circuit NM. The local noise estimate is 

w combined in control circuit CON with the global measure o 2 of the input signal noise to compute by least 
squares regression two coefficients a and b which are respectively applied to the unfiltered and filtered 
signal to generate the output signal z = a*x + b*y- This can be done by multiplying signal x in a first 
multiplier MUL1 by factor a and by multiplying signal y in a second multiplier MUL2 by factor b. The 
outputs of MUL1 and MUL2 are added in adder ADDO, resulting in a noise-reduced output signal z. 

75 Because the kind of filter MF influences the resulting noise reduction and picture quality, advanta- 

geously different filter types can operate in parallel. Then the problem is, how to select the output of the 
best-fitting filter. In Fig. 4 three filters DMF, MED and AV of different kind are operating, whereby weighted 
averages of each two consecutive filter outputs are calculated. 

Each average is obtained as in Fig. 1 by minimising the squared error * 2 between the output signal z and 
20 the noise-free signal u. A respective filtered output signal z lt i = 0 to 2, is selected by examining the 
respective errors: the selected branch corresponds then to the minimum error. This method allows 
switching between filters depending on the statistical properties of the local filtered and unfiltered pictures. 
The three restoration filters are: 

DMF - directional median filter; 
25 MED - median filter; 
AV - average. 

The sliding window 30 depicted in Fig. 3 has a size of 5 pixels by 3 lines. As the processing is not on a 
frame basis but on a field (F1 ; F2) basis, this size corresponds to a region of 5 by 5 pixels for the interlaced 
picture (lines of field F1). 
30 The directional median filter DMF works as follows: 

Let xji, j = 0 to 2, i = 0 to 4, be the pixels taken into account by window 30. The tested filter directions in Fig. 
2 are called & u i =0 to 5. The related gradients are respectively estimated by: 

go = |xoo-xi2| + |x 2 4-xi 2 |; 

35 gi = |Xoi-Xi 2 | + |x 2 3-Xi2|; 

g 2 = |xo 2 -xi 2 | + |x22-xi 2 j; 

g3 = |X03-X1 2 | + |X21-X12|; 

g* = |xo4-xi 2 | + |x 2 o-xi 2 |; 

gs = 1/2*(|Xio-X l 2| + |Xu-Xi 2 | + |Xi3-Xi 2 | + |Xi/v-Xi2|)- 

40 

The selected direction is given by the smallest gradient. The respective output signals yi of filter DMF 
are median values: 

mo = median (xoo, x- t2 , x 2 <0; 
45 mi = median (xoi, x- l2f x 2 3>; 
m 2 = median (xo2, X12, X22); 
m3 = median (xo3t xi 2 , X21 ); 
nru = median (xo* t xi 2 , x 2 o); 
ms = median (xi 1 , xi 2 , X!3>. 

50 

If two or more of these directions result in a minimum gradient, the average of the respective signals is 
taken to generate the output yi . Such multiple minimum gradients appear at less than 1% of pixels. 
The output y 2 of filter MED is: 
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y2 



median (xji) 
0<i<4 
0< j <2 
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The output y3 of filter AV is: 
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y 

0<i<4 
0<j<2 
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75 In Fig. 4 the circuits ADD, DMF, SUB2, CONO, MUL1 , MUL2 and ADDO correspond to the circuits ADD, MF, 
SUB, CON, MUL1, MUL2 and ADD of Fig. 1. The input of filter DMF is also connected to a second filter 
MED and to a third filter AV and can also be connected to an additional filter AF which forms the output 
signal y 0 from the noisy input signal x. The input signal x is also fed to a first, a third and a fifth subtractor 
SUB1, SUB3 and SUBS. The output signal yi , y 2 and y 3 , respectively, of filter DMF. MED and AV, 

20 respectively, is subtracted in the subtractor SUB1 , SUB2 and SUB3, respectively from the input signal x. 
The outputs D 0 . Di , and D2, respectively, of these subtracters are each connected to a respective control 
circuit CONO, CON1 and CON2. The difference signal A 0 . Ai and A 2 . respectively, calculated in a second, 
a fourth and a sixth subtractor SUB2, SUB4 and SUB6, respectively, between signals y 0 ( = x) and yi, yi 
and y 2 . and y 2 and y 3 . respectively, are also fed to the according control circuit CONO. CON1 and CON2. 

25 If filter AF is omitted, SUB1 can also be omitted because Do = An. 

The first control circuit CONO controls the coefficient ao of a first multiplier MUL1 which multiplies signal y 0 
and the coefficient bo of a second multiplier MUL2 which multiplies signal yi. The outputs of these 
multipliers are added in adder ADDO, forming output signal 20. 

The second control circuit CON1 controls the coefficient ai of a third multiplier MUL3 which multiplies 
30 signal yi and the coefficient bi of a fourth multiplier MUL4 which multiplies signal y 2 - The outputs of these 
multipliers are added in adder ADD1, forming output signal . 

The third control circuit CON2 controls the coefficient a 2 of a fifth multiplier MUL5 which multiplies signal y 2 
and the coefficient b? of a sixth multiplier MUL6 which multiplies signal y3. The outputs of these multipliers 
are added in adder ADD2, forming output signal z 2 . 
35 In a seventh, eighth and ninth subtractor SUB7, SUB8 and SUB9, respectively, the difference signals 
between the according output signals zo, zi and z 2 and the input signal x (error on window) are calculated. 
These subtractor output signals are used in a fourth control and mixing circuit COM for connecting one of 
the output signals zo, zi and z 2 with output z, according to the minimum error. 

For calculating the minimum error in the four control circuits, the known noise variance a 2 , as described for 
40 Fig. 1 , is used. Hard or soft switching can be made in control circuit COM. 

The quadratic errors to be minimised in CONO, CON1 and C0N2 are for i = 0 to 2: 



where N R is the total number of pixels in the region R which corresponds to the sliding window 30. The 
notation S u v is used to designate a sum of squares over the region R and is defined as: 



45 




(1) 



50 



= 1/N R * S 
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Su, v 



= £ u(k,l)v(k,l) , 
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where k and I are elements of R and u and v stand either for x, Yj, z t or u. u is not known so that Sy if u 
cannot be computed in such a way. Instead, use is made of an a prior knowledge of the noise statistics. 
It is assumed for simplicity that aj and bj are constant over the region R. Then (equation 1): 

N R s. = a. 2 Z y . 2 + b. 2 E y . +1 2 + 2a. bi l y . y . +i . 2a ^ 

= a i 2 Sy if y ± + b. 2 s yi + 1/y . +l + S M/M ^a.b^y^y.^ ^a.Sy.^ 
-2b.S y . +1 , p 



75 



The constraint is aj + bj = 1. The minimisation with this constraint leads to the following system of 
equations: 

20 aiSyi.yi + biSyi,y j+1 + Xj = $y„u 

aiSyi,y iM + bjSy U i.yi + i + h = Sy iflt u 
B t + bj = 1, 

where \ t is a Lagrange multiplier. 
25 The solution is: 

ai = ( Syj,u -Sy i+1 ,u -Syj.y i+1 + Sy i+1 ,y i+1 )/Aj a, = ( Syi,u -Sy i+1 ,u -Sy i( y j+1 + Sy j+1l y i+1 )/Aj 

bj = (-Sy it u + Sy i+1 ,u -Syj,y j+1 + Sy^ )/A f 

where Aj= Syj.v; -2Syj,y i+1 + Sy j+1 ,y i+1 = Sy r y i+ i,yryi + i 

30 

Because u = x - n, 

Syj.u = Syj.x -Syj,n = Sx.ys -Su f \n -Sn,\n = Sx,Vj - Pi a 2 where for each filter F s : u/ is the filtered noise free 
signal, ny is the filtered noise, o 2 = E(n 2 (k, I)) is the known noise variance and Pi = E(n(k, \)n'(k, I)) is the 
residual noise correlation coefficient after filtering with Fj (i.e. p = 1: no filtering; p = 0: full filtering, no noise 
35 left). 

The assumption of stationary noise allows here easy estimation of these coefficients Pi . But if they are 
known for another kind of noise, it is possible to extend the filter to less restrictive conditions. 
Using these coefficients, a- t and bj can be estimated to: 

40 ^ = ( Sx.Vj -Sx.y i+1 -Syj,y i+1 +Sy j+1 ,y j+1 - Pi a 2 + Pj+1 a 2 )/Aj 
bj = (-Sx,* + Sx,y i+1 -Syi,y i+1 + Sy^j + Pi c 2 - Pi+ ia 2 )/Aj 

or 

45 ai = (Sx-y j+l ,y r y i+1 - a 2 ( Pi - Pj+ i))/Syry i+ 1 ,yry i+ 1 
bj = (Sy r x,y r y j+ i + c 2 ( Pi - PH i))/Sy r y H1 ,y r y i+1 

Therefore a t and bj are known once the correlation term Sx-y i+1 ,yj-y j+l has been computed within the 
control circuits CONO, CON1 and CON2 for the current position of the window 30. The output signals Do, Di 
so and D 2 , respectively, of the according subtractors SUBi, SUB2 and SUB3 are used for computing these 
correlation terms. For good stability of the system it is required that z t is always between yj and y i+1( i.e. 0 
3 a { £ 1 . 

For selecting in control and mixing circuit COM the optimum filter branch, it is possible to derive from 
the values of a and b the quadratic errors: 

55 

N R €j = Sz r u,Zj-u = Szi,Zi -2Szi,u + Sum 
where 
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Su.u. = Sx-n.u. = Sx.u -Sn.u = Sx.x -Sx.n = Sx,x -Sn.n -Su,n = Sx,x -o 2 
and 

5 

Szj.a = Saiyi + biy i+1 ,u = ajSyi.u + bjSy j+1 ,u = ajSx.yj -a<p t Q 2 + bjSx,y i4 ., -bjp i+ io 2 = Sx,Zj -a 2 - 
(a iPi +bip i+ i). 

Therefore 

10 

Np <j = Szi.2j-2Sx.Zj +2a 2 (a t pi +b iPj+1 ) + Sx,x -o 2 = Sz r x,z r x + o 2 (2a jPj + 2b lPi+l -1). 

The correction term e t = o 2 (2ajpj + 2bjp j+ r1) can be explained by the fact that the more a given filter 
smooths the picture, the higher the term Sz r x,z r x will be. Therefore, before the terms generated by filters 
75 of different efficiencies can be compared, they must be compensated by a specific kind of function. 

The graph of this function is given in Fig. 5 to illustrate this effect when using the current filter type (AV, 
MED, DMF, AF) parameters. The correction term ej is shown in relation to coefficient a;. 

The filter AV of type 'average' is well suited for plain areas. The two median filters DMF and MED 
preserve the edges and a good resolution. 
20 To improve the contrast, enhancement processing can be implemented after the noise reduction process- 
ing. 

If the window 30 is enlarged, also other picture degradations can be reduced, for example blurs and stripes. 
Advantageously due to field processing no blurs are introduced on moving objects. 

In order to reduce the complexity and the costs of the hardware implementation it is possible to modify 
25 the described filters in the following way: 

* the size of the window for error estimation can be reduced to 5 pixels by 1 line: therefore no 
additional line memories will be needed (only 2 line memories for the input filters); 

* filter DMF can be a directional median filter with only four directions instead of six: do, d2, d* and ds; 

* filter MED can be a median with a window of three pixels by three lines; 

30 * limitation of the number of bits for the computation of the weighting coefficients (4 bits, for example). 
These simplifications have the following influences on the processed pictures: 

* if the size of the estimation window is smaller, the simplified filter will tend to filter more than the non 
simplified one, particularly when there are local horizontal line structures, because no care is taken of 
the previous and following lines even if they are very different in their local structure; 

35 * simplifications of the input filters do not create too many degradations on the output signal since the 
filtering is adaptive and always the best solution between the three outputs zo, zi and z 2 is chosen; 

* the reduction of the number of bits of the weighting coefficients creates no significant loss of 
accuracy: 1 or 2 units for an amplitude range of 256. 

The big advantage of this new filter is to allow with very few memories and simple hardware implementation 
40 a noise filtering that is adaptive simultaneously locally and globally. 

Because of this double ability of adaptation the filtering is optimised even with simple filters at the input. 

The inventive noise reduction can be used in any video or audio processing devices, e.g. TV receivers 
and VCR's. 

45 Claims 

1. Method for noise reduction, in which for estimating a local noise value within a sliding window (30) of a 
picture the quadratic error between an unfiltered noisy input signal (x) and the filtered (MF) input signal 
(y) is computed and used for calculating (CON) a weighted (a, b) average of the filtered and unfiltered 
50 input signal as output signal (z), characterised in that for calculating (CON; CON0, CON1 , CON2) said 
weighted (a, b; ao, bo, ai, bi, aa, te) average also a global noise value (NM. o 2 ) is used, whereby for 
computing said local noise value pixel difference signals (A; Ao, Ai A2) of said window are calculated 
(CON; CON0, CON1, CON2). 

55 2. Method according to claim 1, characterised in that said filtering is carried out in two or more 
branches of different filter characteristics (DMF, MED, AV) each having a weighted average output 
signal (zo, zi , z 2 ) derived from the filtered input signals (y 0 , yi ; yi , y 2 ; y 2 . y 3 ) of each two consecutive 
branches, whereby for computing said local noise value in each of said branches additionally respective 
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pixel difference signals between each two consecutive of said filtered input signals (yi, V2; y2, yo) are 
used and for obtaining the final output signal (z) that branch output signal (zo. Zi, z 2 ) is selected (COM) 
which has a minimum error compared to said input signal (x). 

5 3. Method according to claim 2. characterised in that when comparing said minimum error, a correction 
term related to the respective filter characteristic is added to each of the branch errors. 

4. Method according to any of claims 1 to 3, characterised in that the sum of the two weighting factors 
(ao, bo; ai t bi ; a 2l D2 ) for said weighted average/s is 1. 

70 

5. Method according to claim 1 and 4, characterised in that the one weighting factor (a) related to said 
unfiltered input signal (x) of the two (a, b) weighting factors for said weighted average is in the range 0 
to 1. 

75 6. Method according to claims 4 and 5 and to claims 2 or 3, characterised in that the respective one (ao, 
ai, a 2 ) of each two (ao, bo; ai, bi ; a 2 , b 2 ) weighting factors for said weighted average is in the range 0 
to 1. 

7. Method according to any of claims 2 to 4 and 6, characterised in that at minimum two of the following 
20 filter characteristics are used in said branches: directional median filter (DMF), median filter (MED), 

average (AV), whereby for said directional median filter at minimum four directions (d 0 -d 5 ) are 
evaluated in said sliding window (30). 

8. Method according to any of claims 1 to 7, characterised in that said sliding window (30) has a size of 
25 five pixels by one to five lines when related to a frame basis. 

9. Apparatus for a method according to any of claims 2 to 4 and 6 to 8, including: 

- a first filter (DMF), a second filter (MED) and a third filter (AV) having a common input signal (x) 
and a different kind of characteristic; 

30 - a first control circuit (CON0) controlling the weighting values (ao, bo) of a first multiplier (MUL1) 

operating on said input signal (x) and a second multiplier (MUL2) operating on the output signal 
(yi ) of said first filter, whereby the outputs of these multipliers are combined in first adding means 
(ADDO); 

- a second control circuit (C0N1) controlling the weighting values (ai , bi) of a third multiplier 
35 (MUL3) operating on the output signal (yi ) of said first filter and a fourth multiplier (MUL4) 

operating on the output signal (y 2 ) of said second filter, whereby the outputs of these multipliers 
are combined in second adding means (ADD1); 

- a third control circuit (CON2) controlling the weighting values (a 2 , b 2 ) of a fifth multiplier (MUL5) 
operating on the output signal (y 2 ) of said second filter and a sixth multiplier (MUL6) operating on 

40 the output signal (y 3 ) of said third filter, whereby the outputs of these multipliers are combined in 

third adding means (ADD2); 

- first subtraction means (SUB2) which feed the difference signal (Ao) between said input signal (x) 
and the output signal (yi ) of said first filter to said first control circuit (CON0); 

- second subtraction means (SUB3) which feed the difference signal (Di ) between said input signal 
45 <x) and the output signal (y 2 ) of said second filter to said second control circuit (CON1); 

- third subtraction means (SUB4) which feed the difference signal (Ai) between the output signal 
(yi) of said first filter and the output signal (y 2 ) of said second filter to said second control circuit 
(CON1); 

- fourth subtraction means (SUBS) which feed the difference signal (D 2 ) between said input signal 
so (x) and the output signal (y 3 ) of said third filter to said third control circuit (CON2); 

- fifth subtraction means (SUB6) which feed the difference signal (A 2 ) between the output signal 
(y 2 ) of said second filter and the output signal (y 3 ) of said third filter to said third control circuit 
(CON2); 

- further subtraction means (SUB7, SUB8, SUB9) which feed the respective difference signals 
55 between said input signal (x) and the output signals of said adding means (ADDO, ADD1 , ADD2) 

to a control and mixing circuit (COM) which selects an output signal (zo, zi, z 2 ) of said adding 
means according to the minimum error derived from these difference signals and according to a 
global noise variance value (a 2 ), 
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whereby for each of the branches said weighting values are derived from error signals calculated in the 
respective control circuit from said difference signals and from said global noise variance value (a 2 ). 

10. Apparatus according to claim 9, characterised in that said first filter (DMF) is a directional median 
5 filter wherein at minimum four directions (do-ds) are evaluated and said second filter (MED) is a median 

filter and said third filter (AV) is an averaging filter. 
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